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1. Introduction 

Recent relativistic heavy ion collision experiments have revealed various remarkable prop- 
erties of QCD at finite temperatures and densities, suggesting the realization of the QCD phase 
transition from the hadronic matter to the quark-gluon plasma (QGP). In order to extract an un- 
ambiguous signal for the transition from the heavy ion collision experiments, it is indispensable 
to make quantitative calculation of the thermal properties of QGP from first principles. Currently, 
the lattice QCD simulation is the only systematic method to do so. By now, most of the lattice 
QCD studies at finite temperature and chemical potential have been performed using staggered- 
type quarks which require less computational costs than others. However, the lattice artifacts of the 
staggered-type quarks are not fully understood. Therefore, it is important to compare the results 
from other lattice quarks to control and estimate the lattice discretization errors. 

We have started systematic studies of finite temperature/density QCD using Iwasaki's RG- 
improved gauge action [jl|] and a clover-improved Wilson quark action In particular, we per- 
form simulations along the lines of constant physics (LCP's) to clearly extract the temperature- and 
density-dependences. As the first project in this direction, we are carrying out simulations of Nf = 2 
QCD on an Nf x N t = 16 3 x 4 lattice at m PS /m Y = 0.65 and 0.80 in the range T /T pc ~ 0.76-3.2, 
where T pc is the pseudocritical point along the LCP Basic properties of the system at finite tem- 
peratures, such as the phase structure and the equation of state, have been studied in |g[ Qj. In 
contrast to the studies with the staggered-type quarks, the expected 0(4) scaling for Nf = 2 QCD 
was confirmed. We extend the study to analyse detailed properties of QGP at finite temperature 
and chemical potential. 

In this paper, we report on the status of our simulations and present a preliminary result for 
the spatial string tension at finite temperature. Results for the free energies of static quarks at zero 
chemical potential are presented in [g|. Preliminary results of the equation of state and susceptibil- 
ities at non-zero chemical potential by the Taylor expansion method are given in [||]. 



2. Lattice action 

We employ the RG-improved gauge action [|l|] and the Nf = 2 clover-improved Wilson quark 
action ^ defined by 

S = S g + S q , (2.1) 
S g = ~PUco t W^W + ci I W;f(*)Y (2.2) 

x \ ^<v;,u,v=l ^v;ii,v=l J 

S q = £ L'/^.v'A'. (2-3) 

f=l,2x,y 

where j3 = 6/g 2 , c\ = —0.331, cq = 1 — 8ci and 

D X y = 8^-KY,{(l-y»)U x ^8 x+flJ + (\ + Y»)Ul^8 x , y+fl }-8 X) .c S wK £ (7 MV F MV . (2.4) 

Here K is the hopping parameter, F^ v is the lattice field strength with f^ v the standard clover- 
shaped combination of gauge links, F^ v = l/Si(f^ v — f^v)- F° r tne clover coefficient csw, we 
adopt a mean field value using W lxl which was calculated in the one-loop perturbation theory [II]], 
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Figure 1 : Phase diagram for RG improved gauge action and clover improved Wilson quark action for N t =4-. 



The phase diagram of this action in the (j8 , K) plane has been obtained by the CP-PACS Collabo- 
ration []3|, |j] as shown in Figjl| The solid line K C (T = 0) with filled circles is the chiral limit where 
pseudoscalar mass vanishes at zero temperature. Above the K C {T = 0) line, the parity-flavor sym- 
metry is spontaneously broken. At finite temperatures, the cusp of the parity-broken phase retracts 
from the large j8 limit to a finite j8. The solid line K C (T > 0) connecting open symbols represents 
the boundary of the parity-broken phase. The red line K t represents the finite temperature pseudo- 
critical line determined from the peak of Polyakov loop susceptibility. This line separates the hot 
phase (the quark-gluon plasma phase) and the cold phase (the hadron phase). The crossing point 
of the K, and the K C (T = 0) lines is the chiral phase transition point. 

3. Determination of lines of constant physics and simulation parameters 

For phenomenological applications, we need to investigate the temperature dependence of 
thermodynamic observables on a line of constant physics (LCP), which we determine by the ratio 
raps / ray of pseudoscalar and vector meson masses. For our purpose, we need LCP in a wider range 
of parameters than [Q]. Therefore, we re-analyze the data for rapsa and my a at zero temperature 
shown in Fig.|| determined by the CP-PACS Collaboration []3| |], 0, ||]]. The colored solid lines in 
Fig.|| shows our results for LCP corresponding to raps/ray = 0.65,0.70,0.75,0.80,0.85,0.90 and 
0.95. The green line denoted as K c represents the critical line, i.e. raps/my = 0. Our LCP's are 
consistent with those of [Q] in the range of the previous study. 

We also reanalyze the lines of constant temperature T /T pc . The temperature T is estimated by 
the zero-temperature vector meson mass mya([},K) using 



csw = (W 



lx i)-3/4 = (i_o.8412 j 8- 1 )- 3 / 4 . 



(2.5) 




1 



\ x mya(p,K) ' 



(3.1) 
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Figure 2: Pseudoscalar meson mass squared (left) and vector meson mass (right) as a function of l/K for 
several values of j3 at T = 0. 




Figure 3: Lines of constant mps/wy an d constant T /T ps in the (j3,K) plane. 



The lines of constant T /T pc is determined as the ratio of T/my to T pc /m\ where T pc /my is ob- 
tained by T ' /my at K t on the same LCP. The bold red line denoted as K t (N t = 4) in Fig.|| rep- 
resents the pseudocritical line T /T pc = 1. The dashed lines represent the results for T/T pc = 
0.8,1.2,1.4,1.6,1.8,2.0 atN t =4. 

We perform finite temperature simulations on a lattice with a temporal extent N t = 4 and a 
spatial extent N s = 16 along the LCP's for mps/my = 0.65 and 0.80. The standard hybrid Monte 
Carlo algorithm is employed to generate full QCD configurations with two flavors of dynamical 
quarks. The length of one trajectory is unity and the molecular dynamics step size is tuned to 
achieve an acceptance rate greater than about 70%. Runs are carried out in the range j8 = 1.50— 
2.30 at twelve values of T/T pc ~ 0.82 — 3.2 for mps/wiy = 0.65 and eleven values of T /T pc ~ 
0.76 — 2.5 for raps/my = 0.80. Our simulation parameters and the corresponding temperatures are 
summarized in Table 1. The number of trajectories for each run after thermalization is 5000 — 6000. 
We measure physical quantities at every 10 trajectories. 
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Table 1: Simulation parameters for m P s/m v = 0.80 (left) and m P s/my = 0.65 (right). 



p 


K 


T/Tpc 


Traj. 


1.50 


0.143480 


0.76446 


5500 


1.60 


0.143749 


0.79544 


6000 


1.70 


0.142871 


0.84346 


6000 


1.80 


0.141139 


0.92507 


6000 


1.85 


0.140070 


0.98642 


6000 


1.90 


0.138817 


1.07619 


6000 


1.95 


0.137716 


1.19836 


6000 


2.00 


0.136931 


1.34778 


5000 


2.10 


0.135860 


1.69025 


5000 


2.20 


0.135010 


2.07325 


5000 


2.30 


0.134194 


2.51093 


5000 



P 


K 


T/T pc 


Traj. 


1.50 


0.150290 


0.82434 


5000 


1.60 


0.150030 


0.86471 


5000 


1.70 


0.148086 


0.94442 


5000 


1.75 


0.146763 


1.00024 


5000 


1.80 


0.145127 


1.07466 


5000 


1.85 


0.143502 


1.17857 


5000 


1.90 


0.141849 


1.31675 


5000 


1.95 


0.140472 


1.48262 


5000 


2.00 


0.139411 


1.66828 


5000 


2.10 


0.137833 


2.09054 


5000 


2.20 


0.136596 


2.59279 


5000 



2.30 0.135492 3.21536 5000 



4. Spatial Wilson Loop 

Using the stored configurations, we are carrying out a series of studies on the nature of the 
quark-gluon plasma. In this report, we present our preliminary results on the confinement in the 
spatial directions at high temperature. 



In previous quenched studies [|9j, |lfj, |llj, [12J], Wilson loops in spatial directions are found to 
show non-vanishing spatial string tension a s even at T > T pc , which is called the spatial confine- 
ment. We study this phenomenon when there exist dynamical quarks in the system. Altough quarks 
are expected to decouple from the spatial observables for T 3> T pc due to dimensional reduction 
and thus do not affect o s in the high temperature limit, it is not obvious whether the same is true 
near T pc . 

As a first trial, we evaluate o s assuming the simplest ansatz for the spatial Wilson loops W(I,J) 
with the size I xj: 

-In W (/,/) = o s IJ+o pen {2I + 2J) +C W , (4.1) 

where o s , o per [ and C w are fit parameters. The results for \/~ojT)/T pc are shown in Fig.|] as a 
function of T/T pc . We find that a/ o s (T)/T pc tends to a non-vanishing value below T pc while it 
increases linearly in T/T pc above T pc . Similar behavior was observed in the quenched case too [§]. 

Let us now compare the results with a calculation of ^/a^ by the three-dimensional effective 
theory assuming the dominance of the gauge part. We would expect the following behaviour for 
the spatial string tension, 



^a s (T) =cg 2 (T) T, (4.2) 
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Figure 4: The square root of the spatial string tension over T pc (left) and T (right) as a function of T/T pc for 
mps/m rm y — 0.65,0.80. The dashed line (right) shows the two-loop fit. 



where g 2 (T) is the two-loop temperature dependent running coupling constant in four dimensions, 



s y ' 48?r 2 A 7687T 4 \ AJ K ' 



We carry out a fit to our data shown in Fig.|], regarding c and A free parameters. From the two 
parameter fit, we find 

c = 0.54(6), A/T pc = 0.14(4). (4.4) 



We note that these values are similar to those obtained in a quenched study Qllj \12\\: c = 0.566(13) 
and A/T pc = 0.104(9). 

We also note that there is a parameter-free three-dimensional effective theory prediction for 
G s (T I T pc , T pc / Aj^g-, Jl/T) [f3|], where Jl is the MS scheme scale parameter. The J[/T is fixed as 
in [13]. We vary the value of T pc /A^ in the range 0.53 — 0.77 [g[ [T^]. The results for y/o^/T are 
compared with our data in Figj5[ While the 1-loop prediction of the effective theory (shown by the 
region bounded by two green lines in Figj5]) deviates from the lattice data, the 2-loop result (shown 
by the region bounded by two black lines) is roughly consistent with our data. A more detailed test 
is left for a future work. 



5. Conclusions 

Most of the lattice QCD studies at finite temperatures and densities have been done using 
staggered-type quarks. To control the lattice artifacts, comparisons with other lattice quarks are 
indispensable. Therefore, we have started a systematic study of QCD at finite temperature and 
density with an improved Wilson quark action. 

As a first step, we performed simulations of Nf = 2 QCD on an A^ 3 x N t = 16 3 x 4 lattice. 
We have identified the lines of constant physics and studied the temperature-dependence of various 
quantities at nips/my = 0.65 and 0.80 in the range T /T pc ~ 0.76-3.2. A preliminary result of the 
spatial string tension o s (T) in the quark-gluon plasma shows a behavior consistent with \J <J S (T) = 
c g 2 (T) T where c takes a value close to that in the quenched case (gluon plasma). o s (T) is also 
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Figure 5: Comparison of our data and 1 -loop and 2-loop predictions from three-dimensional effective theory. 
Jl/T is fixed as in rffjl]. T pc /A^ is varied in the range 0.53 — 0.77. 



consistent with the parameter-free prediction of the three-dimensional effective theory in the 2- 
loop order. Further results on the static quark free energies at finite temperatures for different color 
channels are given in [Q]. Results on the equation of state and susceptibilities at non-zero chemical 
potential with Wilson-type quarks are shown in [0]. 
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